Current distribution inside Py/Cu lateral spin-valve device 
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We have investigated experimentally the non-local voltage signal (NLVS) in the lateral permalloy 
(Py)/Cu/Py spin valve devices with different width of Cu stripes. We found that NLVS strongly 
depends on the distribution of the spin-polarized current inside Cu strip in the vicinity of the 
Py-detector. To explain these data we have developed a diffusion model describing spatial (3D) 
distribution of the spin-polarized current in the device. The results of our calculations show that 
NLVS is decreased by factor of 10 due to spin flip-scattering occurring at Py/Cu interface. The 
interface resistivity on Py/Cu interface is also present, but its contribution to reduction of NLVS is 
minor. We also found that most of the spin-polarized current is injected within the region 30 nm 
from Py-injector/Cu interface. In the area at Py-detector/Cu interface, the spin-polarized current is 
found to flow mainly close on the injector side, with 1/e exponential decay in the magnitude within 
the distance 80 nm. 



PACS Numbers: 75. 70. Pa, 75.70.Kw, 85.70.Kh 



I. INTRODUCTION 

Spintronics is a quickly evolving field providing the 
possibility to manipulate spin degrees of freedom in the 
solid state systems 1,2 . Spin injection, transport and de- 
tection in metals and semiconductors are of particular 
importance to construct effective spintronic devices such 
as a spin battery 3 and spin torque transistor 4 etc. Such 
devices have great advantages over the conventional elec- 
tronic devices because of additional spin functionalities. 
To realize the device it is a key to obtain both large 
spin-polarized current and spin accumulation. It is also 
important to understand the diffusive nature of the spin- 
polarized currents in multiterminal devices. 

Recently, the non-local probing method was proposed 
by Jedema et al. 5 ~ 7 to extract only spin-polarized current 
contribution from the spin-dependent phenomena and to 
reduce spurious effects such as Hall effect and anisotropic 
magnetorcsistancc. They succeeded in detecting the clear 
spin-accumulation signal in the vicinity of the non(N)- 
/ferro(F)-magnetic planar junction by the non-local spin- 
voltage (NLVS) even at room temperature 6,7 . Further- 
more, non-local technique maybe useful to induce spin- 
injection magnetization reversal 8 without the charge cur- 
rent, leading to the solution for the energy dissipation 
problem due to Joule heat. 

Here, we study experimentally the distribution of the 
spin-polarized current in non-local configuration. So far, 
the spin-polarized current transport is analytically inves- 
tigated using one dimensional (ID) Boltzmann diffusion 
model 6,9-12 . As these models predict too large NLVS, we 
have developed formalism to calculate spatial (3D) distri- 
bution of spin-polarized current. However, large decrease 
of NLVS can not be attributed to spatial distribution of 



spin-polarized current, and we attribute it 
tcring at Py/Cu interface. 



to spin scat- 



II. DEVICE PREPARATION AND 
EXPERIMENTAL RESULTS 

We fabricated lateral spin-valve devices consisting of 
two Py wires bridged by Cu strip by means of nano- 
fabrication techniques. Figure 1 shows a scanning- 
electron-microscope (SEM) image for one of the fabri- 
cated devices. First, we fabricated both Py wires of width 



w Py 



120 nm and of thickness ip y = 20 nm with the 



spacing of ^cu = 170nm by electron-beam lithography 
and lift-off technique. Py layer was evaporated by an 
electron-beam gun at 2 xl0~ 8 Torr. Ends of the first 
Py wire are connected to large pads pattern for assisting 
the nucleation of the domain wall, although ends of the 
second one are flat-end shaped. Hence, each Py wire has 
different switching field. 

Both Py wires are bridged by Cu strip of thickness 




FIG. 1: SEM image of the fabricated lateral spin-valve device. 
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FIG. 2: DNLVS for system with 
(a) 2-wires and (b) with 3-wires. 
(c) Difference of local voltage sig- 
nal (DLVS) for system with 2 
wires. 
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FIG. 3: (a)(b) SEM image of a de- 
tail of the lateral spin-valve device 
with iajcu = 100, 500 nm, respec- 
500 tively, with sketched current flows 
and "cross" and "half" detection 
configuration, (c-f ) NLVS as a func- 
tion of external magnetic field, ob- 
tained for "cross" and "half" con- 
figuration for ibcu = 100, 500 nm. 
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tc-a = 80 nm, having widths w'cu =100, 300 and 500 nm 
for three different devices. Prior to Cu deposition, the Py 
surface was cleaned by Ar + bombardment and then sam- 
ple was shortly taken out-of-vacuum to change vacuum 
chamber. Then Cu was evaporated by resistance heat- 
ing. The contact resistance of the interface was found 
ohmic and very low indicating a transparent contact. 
The conductivity of Cu is cr C u,RT = 48.1 x 10 6 f2 _1 m _1 , 
ccu,4K = 131 x 10 6 f2 _1 m _1 at room temperature, 4 K, 
respectively. 

Notice that the present Cu strip has smaller residual 
resistivity than that of Jedema at el. 6 . The NLVS mea- 
surements were performed at room temperature with the 
magnetic field applied parallel along the Py-wires. 

NLVS measurements were performed using a standard 
current-bias lock-in technique at room temperature. We 
measured the NLVS as a function of external magnetic 
field by using two different probe configurations, called 
"half" and "cross" . The difference between both config- 
urations is whether the current and voltage probes are 
located on the same side or not as sketched in Figs. 2(a) 
or 3(a). The one-dimensional (ID) diffusion model 12 ' 13 
predicts that the obtained NLVS should be the same. 
However, as the spin-polarized current has the spatial dis- 
tribution, the NLVS shows the difference between both 
probe configurations. 

Figures 3(c) and 3(e) show a NLVS for wc u = 100 nm 
with "half" and "cross" probe configuration, respectively. 
The obtained difference of NLVS between parallel and an- 
tiparallcl magnetization (DNLVS) is 0.7 mfi and 0.6 mfi 



at room temperature, respectively Figures 3(d) and 3(f) 
show NLVS for wc u = 500 nm with "half" and "cross" 
configurations, providing 0.6 and 0.1 mf2, respectively. 

Experimental values of DNLVS as a function of u>cu are 
presented on Fig. 4(a). Experimental data show that the 
difference between "cross' and "half" in DNLVS increases 
with increasing wqu- 

The other parameters (at room temperature) used in 
our calculations are as follow: Py conductivity op y = 
7.3 x 10 6 f2 _1 m _1 , Py bulk spin asymmetry coefficient 
= 0.7 (<7 T)Py = a Py (l + /?)/2, a UPy = a Py (l - 0)/2) 
(Refs. 14 ~ 16 ), spin-flip lengths A Py = 4.3 nm (Ref. 14 ) and 
Acu = 350 nm (Ref. 13 ). Py wires have width iyp y = 
120 nm, thickness tp y = 20 nm and separated by distance 
(cu = 170 nm. The Cu strip has thickness £cu = 80 nm 
with widths wcu = 100, 300, 500 nm. 



III. ID CALCULATIONS OF DNLVS 

In the literature, there are two models describing NLVS 
(and DNLVS) inside metallic lateral spin-valve device: 
one given by Jedema et. al. 13 and the other by Takahashi 
and Maekawa 12 . Both models approximate the device 
into ID wire circuit, in which the propagation of electro- 
chemical potential /i| / j and spin-polarized current Jj n 
is described by standard Valet-Fert model 10 . At an in- 
tersection point of several wires (hereafter called node), 
e.g. intersection of Cu and Py wires, the boundary con- 
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FIG. 4: Experimental value of DNLVS as a function of Cu 
wire width wcu compared with ID and 3D models. Exper- 
imental data (square) compared with ID and 3D models. 
R* and R s are interface and surface scattering resistances 
for Py/Cu interface, units in film 2 . For detail see Sec. Ill 
and IV B. 

ditions, expressed as generalized Kirchoff's laws are 

H J »d/i =° 

/VT/i = const T/1 

where n is an index of all the wires connected to a given 
node. Hence, the /if n (which can be understand as 
a voltage here) is the same for each ID wire attached to a 
given node, and Jy± is conserved while flowing through 
each node. 

The model of Jedcma et al. 13 has two assumptions, 
which arc not fulfilled in our case of Py/Cu device: (a) 
they assume cross-sectional areas of all the wires in the 
device were the same (i.e. they considered more continu- 
ity of up and down current densities j\ /I than up and 
down currents Jj n at each node) and (b) they assume 
(\f,Xn) 3> (wf,wn), where (wf,wn) are widths of F, 
N wires, respectively. The comparison of DNLVS ob- 
tained from this model (when extended to the case for 
different cross-sectional areas of wires) with our exper- 
imental data is shown on Fig. 4(b) (dashed-dot line), 
showing that this model predicts about 40 x large value 
than experimental one. 

These drawbacks were partly overcome by Takahashi 
and Maekawa 12 , assuming that (a) \p <C (wf,wn) <C 
An and (b) that current at F/N interface is homoge- 
neous. Later we will show that assumption (b) is not 




FIG. 5: The sketch of 3D network of spin-dependent- 
resistance-elements (SDRE). Circle inset sketches that SDRE 
consists of spin-up and spin-down resistances and of the shunt- 
ing resistors between up and down channels. Note that each 
node and wire on the sketch represents a "bus" containing 
spin-up and spin-down channels. 

correct for ohmic junctions, but is correct for tunnel 
junctions. Although they derived their model from basic 
equations, the same results can be obtained when both F- 
injector and F-detector, attached to N wire, are described 
by a standard ID model, where F-wires have effective 
cross-section area as of Py/Cu interface, i.e. in our case 
§f = wpWN- The DNLVS calculated from this model 
is presented on Fig. 4(b) for the case with interface re- 
sistance Rp y / Cu = (solid line) and -Rp y / Cu = 0.5ffhii 2 , 

7 = 0.7 (Ref. 15 ) (dashed line). Note that i?p y/Cu = 

2^p y /cu(l " 7), 4 y /cu = 2Ai^ y/Cu (l + 7 )- This ID 
model describes quite well the experimentally observed 
DNLVS but gives about 10 times larger magnitude than 
the experimental results. 



IV. 3D CALCULATION OF SPIN-POLARIZED 
CURRENT AND ELECTROCHEMICAL 
POTENTIAL 

In order to understand the spin-polarized currents in- 
side device in detail, we have developed model calculating 
3D distribution of /i| / j and spin-polarized current density 
jj/l inside the device 17 . Our model is based on the 3D 
electrical network of spin-dependent-resistance-elements 
(SDRE) (Fig. 5). The response of each SDRE is deter- 
mined by ID models 10,13 . As sketched in inset of Fig. 5, 
each SDRE consists of resistance for spin-up, spin-down 
channels and spin-flip resistance shunting up and down 
channels. This shunting resistance can be regarded as the 
"probability" that electron spins arc flipped when passing 
SDRE. A boundary conditions at each node connecting 
SDRE are given by Eq. (1). 

In this model we can also account surface or interface 
resistance (scattering), AR SS or AR Sl respectively, short- 
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FIG. 6: The profile of j sp = jj — jj through 
Cu/Py(20)/Cu(20)/Py(20)/Cu pillar structure, dimensions 
in nm, calculated for ID VF model (full line) and compared 
with our 3D calculations with perpendicular-to-interface grid 
size 1 nm, 5 nm, 10 nm. Lateral (parallel-to-interface) grid 
size is 10 nm. 

cutting up and down channels at the surface or interface. 
For detailed description of the formalism please see 17 . 



A. Plausibility of 3D calculations 



FIG. 7: Dependence of DNLVS on lateral grid size. Details 
in Sec. IV A. 



vicinity of Py/Cu interface lateral grid size 5nm. This 
grid configuration is denoted by a circle on the Fig. 7, 
providing agreement with DNLVS converged value 13%, 
in agreement with above discussion. Unfortunately, in 
our calculations, we can not use smaller grid size due to 
numerical limitations. We conclude (i) precision of our 
3D calculations is about 20% (ii) calculated DNLVS has 
tendancy to be overestimated. 



To estimate precision of our 3D calculations, we have 
calculated jyi, Hyi and magnetoresistivity ratio (MR) 
in Cu/Py(20)/Cu(20)/Py(20)/Cu multilayer structure 
(dimensions in nm) using different grid sizes for Py. The 
results of these calculations should be identical with ID 
Valet-Fert model 10 . We investigate calculation precision 
only with grid size of Py, as Ap y <C Xq u ■ 

Figure 6 shows profile of spin-polarized current j sp = 
it — if through antiparallel Py/Cu multilayer structure. 
Used lateral grid size (i.e. grid distances parallel with 
Py/Cu interfaces) is 10 nm, perpendicular grid (i.e. grid 
perpendicular to interfaces) is 1, 5, 10 nm, giving j sp pre- 
cision inside Py being 4%, 9%, 15%, respectively. 3D 
calculations gives larger value of MR by 8%, 16%, 33% 
than ID calculation. For in-plane grid size 5 nm, the MR 
is larger by 4%, 11%, 26%. It shows that (i) with de- 
crease of grid size, j sp and MR converge to correct values 
and (ii) small perpendicular grid size is more important 
than in-plane one. 

Figure 7 shows a dependence of DNLVS in the lateral 
spin-valve structure on various lateral (i.e. parallel with 
substrate surface) grid size. The simulated device is dif- 
ferent than real one; two Py wires of 15 nm-thick and 
50 nm- width are separated by a distance of 80 nm and 
bridged by 55 nm-thick, 50 nm-width Cu strip. 

DNLVS has been calculated for perpendicular grid size 
5nm (square in Fig. 7) and 2.5 nm (diamond), providing 
larger DNLVS by 18%, 11%, respectively, with respect to 
the converged DNLVS value. In both cases, larger grid 
size leads to larger DNLVS. 

In all simulations of real structure, we used perpen- 
dicular grid size 5 nm, lateral grid size 10 nm and in the 



B. 3D calculations of DNLVS 

Figure 4(b) presents DNLVS calculated from 3D model for 
ARp y/Cu = (circle) and for ARp y/Cu = 0.5ffmi 2 , 7 = 0.7 
(Ref. 15 ) (diamond). Both DNLVS have about the same shape 
and slightly smaller magnitude compared to values from ID 
model by Takahashi and Maekawa 12 (solid and dashed lines 
without symbols). In agreement with experiment, for larger 
wcu, the DNLVS has different values in half and cross con- 
figuration, reflecting inhomogeneous j sp at position of the de- 
tector. As will be shown in Sec. V, the j sp is also strongly 
inhomogeneous at injector position. However, the approxi- 
mative agreement between 3D models and ID models shows, 
that influence of inhomogeneous current injection is not very 
important to magnitude of DNLVS. As our 3D models tends 
to overestimate DNLVS, we conclude that non-local current 
injection decrease DNLVS, but only about 20%. 

Both ID and 3D models shows that the presence of -R Py / Cu 
together with large positive value of 7 increases DNLVS. 
When i? Py/ , Cu > 0, 7 = 0, DNLVS decreases. It may be possi- 
ble that 7 > 0, but this contribution to DNLVS is smeared by 
other contribution decreasing DNLVS. Therefore, in following 
we assume 7 = 0. 

Now, let us discuss which mechanism decrease DNLVS. To 
be more sure with analysis, we take into account more experi- 
mental data (which are going to be published elsewhere 18 ) on 
two different samples, fabricated exactly by a way as previous 
sample. 

• 3- wires system consisting of two Py wires of 20 nm- 
thick, 100 nm-width, separated by a distance of 400 nm 
and bridged by 80 nm-thick 100 nm-width Cu strip. Be- 
tween both Py wires, there is third 100 nm-width wire 
[Fig. 2(b)], consisting either of Cu (having thickness 
80 nm), or Py wire (having thickness 20 nm), or there 
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0.7 
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0.77 
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0.78 
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DNLVScross, WCu = 100 
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0.76 
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DNLVShaif, ™ C u = 300 


0.6 
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0.87 
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0.40 


0.29 


DNLVScross, WCu = 300 


0.3 


1.2 


0.07 


0.24 


0.08 


0.11 


0.14 


0.21 


0.25 


3- wires, no middle 


0.25 


3.72 


0.034 


0.032 


0.45 


0.47 


0.46 


0.54 


0.48 


3- wires, Cu middle 


0.18 


2.53 


0.025 


0.033 


0.32 


0.33 


0.31 


0.35 


0.28 


3- wires, Py middle 


0.04 


0.76 


0.019 


0.018 


0.037 


0.052 


0.075 


0.16 


0.33 


DNLVShaif, wcu = 250 
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0.51 


2.0 


0.61 


0.29 


0.29 


0.25 


0.27 


0.25 


DNLVS cr0 ss, wcu = 250 


0.4 


1.3 


0.20 


0.33 


0.12 


0.14 


0.15 


0.21 


0.24 


DLVS 2Py , w Cu = 250 


1 


3.8 


2.0 


2.2 


0.50 


0.50 


0.46 


0.51 


0.50 


DLVSiPy, wcu = 250 


0.4 


2.4 


1.7 


1.8 


0.34 


0.32 


0.28 


0.29 


0.26 



TABLE I: Experimental values of DNLVS (in mfi) for various sample structures compared with values calculated by 3D models 
taking into account different processes decreasing DNLVS. Units of resistances are in fOm 2 , units of wire widths in nm. For 
details see Sec. IV B. 



is no third- wire. 

• system consisting of two Py 20 nm-thick wires with dif- 
ferent widths (200nm-width of injector and 100 nm- 
width of detector), separated by 200 nm and bridged 
by 250nm-width and 80 nm-thick Cu strip. In this de- 
vice, we measured both DNLVS and difference of local 
voltage signal between parallel and antiparallel state 
(DLVS) [Fig. 2(c)]. In DLVS case, charge current flows 
through both Py wires. 

In the following, we will discuss possible contributions com- 
ing from (i) surface scattering on Cu, j4_R S s,Cu (ii) surface 
scattering at Py/Cu interface AR StPy /cu (hi) interface non- 
polarized resistance, AR Py ^ Cu . The possible magnitude of 
each contribution has been determined to fit DNLVS for 
wcu = 100 nm and then compared with other experimental 
data. All experimental data and calculated values are sum- 
marized in Tab. I. 



1. Surface scattering on Cu 

Surface scattering on Cu is introduced by an resistance 
^4^ss,Cu shortcutting up and down channel on the Cu surface. 
To decrease DNLVS for ioc u = 100 nm to experimental value, 
0.7mf2, Cu surface scattering has to be -R ss ,Cu = 0.15 fffrn 2 
(when surface scattering is assumed on both side sides and top 
and bottom surface of Cu wire) or -R S s,Cu,sidc = 0.065 ffim 2 
(when surface scattering is assumed to be only on both sides of 



Cu). However, using those surface scattering resistances, the 
DNLVS calculated for 3-wires system (Tab. I) are too small 
compared with experiment, showing that this contribution is 
not a dominant one. 



2. Interface scattering on Py/Cu interface 

The properties of Py/Cu interface is here described by a 
presence of the interface layer, which has its own thickness 
ti, spin- flip- length Aj and conductivity o~i, spin-polarization 
7/ 19 . However, the interface properties should not depend 
on ti (this value is given ad-hoc and is assumed as 1 nm in 
our calculations). Therefore, it is profitable to express inter- 
face properties by Si = ti/Xi and AR Py / Cu = ti/ai, which 
are independent on t/ 19 . Physical meaning of AR* is clear: 
2i?*(l — 7/), 2R*(1 + ji) is a resistance of channel up, down 
through interface layer, respectively. As physical meaning of 
Si is not so clear, we prefer to describe spin-flip scattering by 



interface scattering resistivity 



AR S = AR* 



S sinh S ' 



(2) 



which means a resistance shortcutting up and down channels 
on the interface. 

To decrease DNLVS to experimental value at wcu = 
100 nm, different pairs of AR Py / Cu , ^4-R s ,Py/Cu can be used, 
as shown in Table I. When there is no interface resistance 
(AR Py/Cu = 0), then AR SjPy / Cu = 2.6fffrn 2 . On the other 
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hand, when ^4-R s ,p y /cu = hif then AR Py / Cu — 15 film 2 . Both 
^4-Rs,p y /Cu and ARp y / Cu contribute to decrease of DNLVS. 

Table I and Figure 4(a) shows that none combination of 
pairs AR Py / Cu , -R s ,p y /Cu describes perfectly all experimental 
values, however the agreement with all experimental data is 
within factor of 2-3. Figure 4(a) shows that with increasing 
value of AR Py / Cu , the difference between "half" and "cross" 
DNLVS is reducing, reflecting more homogeneous injection of 
j sp over Py-inj/Cu interface. 

The most relevant interface properties is a pair of values 
AR Py/Cu = lfHrn 2 , i^,p y/ cu = 3.8 fOm 2 (<5 Py /cu = 0.95) 
as for this pair the mutual ratio between DNLVS's for 3-wire 
system (when middle wire is Cu, Py and nothing) agrees with 
experiment. Then all calculated values for 3-wires systems are 
about 1.8x larger then experimental one. The disagreement 
by factor 1.8 x can be related to smaller value of Acu than ex- 
pected 350 nm. The 3-wire configuration with middle Py wire 
is particularly sensitive to AR Py ^ Cll , as its value determines, 
how large amount of j sp is absorbed by the middle Py wire. 

Table I shows that experimental value of DNLVS at racu = 
300 nm is larger than calculated one, particularly for cross 
configuration (experimental DNLVScross =0.3 mil, but calcu- 
lated 0.14 mf2). In another words, DNLVS(uicu) decreases 
slower for experiment than for calculated value. It is prob- 
ably due to presence of a charge current j c h at a position 
of Py-detector for wider iucu, as will be shown in Sec. VB. 
Non-zero j c h inside detector probably causes some additive 
contribution to DNLVS, either due to AMR, either due to the 
scattering related with currents in-plane (CIP), i.e. currents 
flowing parallel with Py/Cu interface 

Figure 4(b) also contain a dependence DNLVS(wcu) cal- 
culated from extended model of Takahashi and Maekawa for 
AR Py/Cu = lfilm 2 , -R s , Py /cu = 3.8fffrn 2 . We can see that 
there is a good agreement with 3D calculations. It shows when 
jsp is homogeneous on detector position, this model predicts 
a correct value of DNLVS. 

The last part of Table I shows an agreement between ex- 
perimental and calculated values of DNLVS and DLVS, de- 
termined for ibcu = 250 nm. We can see that for ARp y / CvL — 
lfffrn 2 , -R s , p y /Cu = 3.8 film 2 , all experimental values are 
about twice large compared to calculated one. Probably, 
here play role similar effects as discussed for DNLVS for 
wcu = 300 nm, as in this case calculated DNLVS is also twice 
smaller than experimental one. 

Resistance AR Py / Cu = lK7m 2 is equal to resistance of 
48 nm of Cu or 7.3 nm of Py. Furthermore, interface scat- 
tering ^4-R 3 ,p y /cu = 3.8H7m 2 corresponds to scattering by Cu 
at length 950 nm and at length 2.5 nm inside Py [Eq. (2)]. Es- 
pecially second value shows that interface scattering is not so 
large, however, it is enough to decrease DNLVS by one order 
of magnitude. 

In conclusion of this Section, we have shown that major 
contribution to small DNLVS is due to interface scattering 
resistance AR StPy /c u , shortcutting up and down channels at 
Py/Cu interfaces. The interface resistivity AR Py / CvL is also 
presented, but its contribution to decrease of DNLVS is only 
minor one. Such a large interface spin-scattering has not been 
observed in 14,15 . It can be related with two factors: 

1. Quality of our Py/Cu interface is lower than in 14,15 . In 
our fabrication process, there are two steps which could 
decrease interface quality. On top of Py we deposited 
and removed photoresist to pattern Cu wire. Before Cu 
deposition, the surface was cleaned by Ar + bombard- 



z» y 




xy cut 



yz cut xz cut 



FIG. 8: The sketch of the device with indicated cut planes. 
The yz, xz' cuts are taken in the center of Py-injector, Cu 
wire, respectively. The xy cut is located 12.5 nm from device 
top, the xz cut is located 7.5 nm from the side of Cu wire. 



ment and then device shortly taken out of vacuum to 
change vacuum chamber. Notice, that Jedema et al. 13 
has used very similar fabrication process as we did. 

2. Contribution of Py/Cu interface spin scattering is miss- 
ing in previous works, investigating Py/Cu system by 
means of magnetoresistivity ratio (MR) 14,15 . Note that 
MR is sensitive to value of j ap passing free layer rather 
that to value of spin accumulation A/i at position of 
free layer 17 . As we have shown 17 , system can provide 
large MR (when large j sp flow through free layer) al- 
though A/i at position of free layer can vanish. When 
A/j, vanishes, then shortcutting of up and down chan- 
nels takes no effect and so interface spin-scattering does 
not occurs at the interface. In such a case, the MR (up 
to some limit) is insensitive to spin-scattering on the 
free-layer / non-magnetic-layer interface. 

On the other hand, non-local technique is particularly 
sensitive to A/j, at detector/non-magnetic- metal inter- 
face. When interface spin-scattering is presented in this 
case, it significantly reduces DNLVS. 

Hence, it may be possible, that small interface scatter- 
ing is presented in both MR and non-local measure- 
ments, but did not take a place in case of MR measure- 
ments. 



V. CURRENT FLOWS INSIDE LATERAL 
SPIN- VALVE STRUCTURE 

In this Section, we present in detail the current inhomo- 
geneity inside lateral spin-valve structure. Figure 8 shows 
a sketch of the device with indicated cut planes, on which 
the calculated current densities are presented on Figures 9, 
11, 12, 13 and discussed in following Sections VA-VC. The 
presented current densities were calculated for parallel mag- 
netizations and for our best interface description AR Py / Cu = 
lffhn 2 , i? Sj p y /cu = 3.8fffrn 2 . For antiparallel magnetiza- 
tions, we get very similar current flows as for parallel one. 
This is in agreement with ID models of non-local devices 12 ' 13 , 
where current flows are exactly the same for parallel and an- 
tiparallel magnetic states. 
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FIG. 9: The yz cut (denned on 
Fig. 8) of the current density of (a) 
current polarized up j] (b) current 
polarized down ]\ (c) charge cur- 
rent j c h = j'f + ji and (d) spin- 
polarized current j sp — jf — jj 
in the device with wcu = 100 nm 
with parallel magnetization. The 
length of arrow is proportional to 
value of a given current; this scal- 
ing is the same for all cuts. 
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FIG. 10: The profiles of j c h (open symbols) and j sp (solid sym- 
bols) taken at the intersection of Py-inj ector/Cu interface and 
the yz cut (Fig. 8). The profiles were calculated for wcu = 
100 nm, for current 1=1 mA, for AR* nt = 0, AR 3 , Py /cu = inf 
(circle) and for _Rf nt = lffim 2 , AR aPy /c u = 3.8 film 2 (tri- 
angles and diamond). j sp and j c h are the same for parallel 
and anti-parallel magnetizations. The symbols denote grid, 
for which the current densities are calculated. 



A. Current description near Py-injector 

Figure 9 shows current density on the yz cut, which is taken 
in the center of the Py-injector wire (Fig. 8). Cuts (a-d) 
correspond to the cases for up and down current densities 
31 1 31 1 respectively, for charge current density j c h = 31 + 31 
and for spin-polarized current density j sp = j j — j\ . All cuts 
show that the current is injected rather sharply through Py- 
injector/Cu interface and then quickly spreads into the whole 
volume of Cu wire. 

The values of j c h and j sp at the intersection of yz cut and 
Py-inj ector/Cu interface are presented on Fig. 10. The profile 
is shown for the device with AR Py ^ Cu = 0, AR aPy / C u = inf 
(circles) and for AR Py / Cu = lK7m 2 , i? s ,p y /cu = 3.8ff2m 2 
(triangles and diamond). Due to spin-flip-scattering on the in- 
terface, the j S p flowing to the interface from Py side (triangles- 
down) is about twice larger than from j sp outgoing the inter- 
face at Cu side (triangles- up). 

It is shown that both j c h (open symbols) and j sp (solid sym- 



bols) are sharply injected within the distance of 25 nm, 35 nm 
from the Py/Cu edge for A7ip y / Cu = lffim 2 , AR a:Py / C u = 
3.8 film 2 and AR Py / Cu = 0, j4-R s ,F y /cu = m f respectively. 
This different 'length-of-injection' is only due to different val- 
ues of AR Py i Cu , and is nearly independent on AR aPy /c u . 
When R Py / Cu is large then obviously the current is more 
spread over the interface and for tunnel contacts is can be 
considered as homogeneous. Furthermore, j ap is positive only 
in the distance of 25 nm or 35 nm from the Py-injector/Cu 
edge, and then its value becomes negative. This means that in 
this region the injector reabsorbs a small part of the injected 
spin-polarized current, which decreases the spin-injection ef- 
ficiency. 

For different values of wcu, the 'length-of-injections' are 
very similar to those presented in Fig. 10. It should be noticed 
that this sharp injection occurs in consequence of small Py 
conductivity, a Py <g trcu and small thickness of Py wire t Py < 
(wcu,tcu)- In other words, larger t Py increases homogeneity 
of the injected current. 



B. Top view on Cu 

Figure 11 presents current density on xy cut (defined in 
Fig. 8) taken at the depth of 12.5 nm from top surface of the 
Cu wire. As already discussed, current is sharply injected at 
Cu/Py-injector edge and hence jf and j sp spread into the Cu 
wire from this edge [Fig. ll(a)(d)]. 

Figure 11(d) also shows that for wcu = 100 nm, j sp at the 
position of the detector is fairly uniform in the y-direction, i.e. 
in the direction parallel to the Py wire. When j'j reaches de- 
tector, it is successively spin-scattered due to very short spin- 
diffusion length Ap y and then current flows homogeneously 
back as j x [Fig. ll(a)(b)]. 

Due to the sharp current injection, j c h makes a whirl in the 
'diffusive' part of the Cu wire, where no charge current was 
expected [Figure 11(c)]. In the present case (wcu — 100 nm), 
the value of j c h originating from this whirl at detector position 
is negligible compared to j sp . However when leu < men (leu 
being distance between Py wires), then j c h ^ j sp at detector 
position. This can be seen on Fig. 12(a) for idcu = 300 nm. 
Fig. 12(b) shows that also j sp for i»c u = 300 nm is inhomoge- 
neous at detector position, having maximal value at one side 
of Cu/Py-detector interface. This explains different values 
between "cross" and "half" configurations. 
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FIG. 11: The xy cut (denned on 
Fig. 8) of (a) i T (b) u (c) j ch 
and (d) j sp in the device with 
wcu = 100 nm, parallel magneti- 
zation. The arrows have the same 
scaling for all cuts, and they are 
magnified on cut's right sides. 
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FIG. 12: The xy cut (defined on Fig. 8) of (a) j ch and (b) 
j sp in the device with uicu = 300 nm. The arrows have the 
same scaling for all cuts, and they are magnified on cut's right 
sides. 



C. Current description near Py-detector 

Figure 13 presents current density on the xz cut (defined 
in Fig. 8) which is taken 7.5 nm from the side of Cu wire. 
Figure 13(a) shows that the flow of ji into detector is also in- 
homogeneous and is dominant at the side of Py-detector/Cu, 
which is close to the injector. As already mentioned above, 
due to very short spin-diffusion length Ap y = 4.3 nm, 
flowing into the detector is immediately reversed inside Py- 
detector and coming back as j± [Fig. 13(b)]. This can be 
understand as a resistance shunting (or "short-cutting") the 
up and down channels. This also explains the behaviour of 
j B p [Fig. 13(d)], whose flow is absorbed by the detector. 

Figure 14 shows and j ap on the intersection between Py- 
detector/Cu interface and xz' cut (defined on Fig. 8), which 
is taken at the center of Cu wire. The vertical dash-dot lines 
show the position of edges inside Py-detector wire embedded 
in the Cu wire, i.e. ranges x £ (—20, 0) and x £ (120, 140) cor- 
responds to side part of Py wire, although range x G (0, 120) 
represents bottom part of Py-detector wire. We can see that 
both jsp and j c h are inhomogeneous, decaying approximately 
exponentially with 1/e decrease length 80 nm. This decay is 
mainly result of the competition between Cu conductivity ucu 
and spin-flip scattering inside Py and Py/Cu interface. 



When the interface resistances AR 



Py/Cu 



= 1 ffW 



AR 



s, Py/Cu 



3.8fffrn are introduced, j Bp flowing to the de- 
tector is decreased (and hence DNLVS is decreased) , as can be 
seen on Fig. 14. Due to presence of ^4-R s ,p y /cu, the current 
flowing to the interface from Cu side (triangle-up) is about 
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FIG. 13: The xz cut (defined on 
Fig. 8) of (a) j T (b) j i (c) j ch 
and (d) j sp in the device with 
aiCu = 100 nm, parallel magneti- 
zation. The arrows have the same 
scaling for all cuts, and they are 
magnified on cut's right sides. 
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FIG. 14: The profile of j c h and j ap on the intersection of Py- 
detector/Cu interface and the xz' cut (defined on Fig. 8) for 
parallel magnetization state. The rest as in Fig. 10. 



3x larger than one outgoing to the injector (triangle down). 
Hence, 2/3 of j sp entering detector are shortcut, absorbed by 
a interface. 

There is also j c h at Py-detector/Cu interface, having value 
about 10% of j S p. The j c h originates because jj and j\ are 
injected to/ejected from the Py-detector at slightly different 
position, i.e. j c h has negative value around x > and positive 
at x < 0. It means that part of jj current, which is injected 
to Py-detector from the side of the Py wire, is ejected as j± 
from its top part. 



VI. CONCLUSION 

We have fabricated lateral spin- valve devices consisting of 
the permalloy (Py) and Cu wires. We have observed that the 
difference of the non-local voltage signal (DNLVS) between 



parallel and antiparallel magnetization has different values 
for "half and "cross" configurations. The difference between 
these two configurations increases when the width of the Cu 
stripe increases. 

To understand observed behaviour in detail, we have de- 
veloped formalism calculating spatial (3D) distribution of the 
spin-polarized current and electrochemical potential inside 
the device. We found that the current distribution inside 
lateral spin-valve device with ohmic-contact is rather com- 
plex interplay between geometry and electrical properties of 
all the involved materials. 

Despite of those large current inhomogeneities, the DNLVS 
calculated from our 3D model are in a good agreement with 
ID model given by Takahashi et Maekawa 12 . However, both 
ID and 3D predicts about 10 x larger DNLVS than experi- 
mental values. We have attributed the smallness of DNLVS 
to interface scattering resistance AR 3 Py /c u = 3.8ffmr 2 short- 
cutting up and down channels at Py/Cu interface. On one 
hand, this value of AR St p y /cn decrease DNLVS by factor of 
10. On the other hand, it corresponds only to scattering 
which occurs inside Py on distance of 2.5 nm. When this in- 
terface scattering resistance can be reduced, DNLVS may be 
enhanced significantly. The fact, that such a interface resis- 
tivity has not been observed before 14,15 may be related either 
to lower quality of our interface, either to insensitivity of MR 
to small surface scattering in some cases. 

Interface resistance ARp y / Cu = lfilm 2 is also presented at 
Py/Cu interface, but its contribution to smallness of DNLVS 
is minor. The value of this resistance mainly modifies the 
j sp inhomogeneity in the structure. Using this description 
of Py/Cu interface, we found agreement with all our experi- 
mental data (local and non-local voltage signals measured on 
systems with two or three Py wires) within factor of two. 

The current is injected from Py-injector to Cu sharply, 
within the distance of 30 nm. Part of the injected spin- 
polarized current is reabsorbed by injector itself. Current flow 
over Py-detector/Cu interface is also inhomogeneous, having 
the largest value on the side of Py-detector close to injector 
and decaying approximately exponentially with 1/e decrease 
within the distance 80 nm. 
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